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The neutrinoless double beta (Ov/3/3) decay process could provide crucial information to determine 
the absolute scale of the neutrino masses, and it is the only one that can establish whether neutrino 
is a Dirac or a Majorana particle. A key ingredient for extracting the absolute neutrino masses from 
0is(3(3 decay experiments is a precise knowledge of the nuclear matrix elements (NME) describing 
the half-life of this process. We developed a new shell model approach for computing the Oz^/3/3 
decay NME, and we used it to analyze the Ov/3/3 mode of 48 Ca. The dependence of the NME on 
the short range correlations parameters, on the average energy of the intermediate states, on the 
finite-size cutoff parameters, and on the effective interaction used for the many-body calculations is 
discussed. 

PACS numbers: 23.40.Bw, 21.60.Cs, 23.40.-s, 14.60.Pq 



I. INTRODUCTION 



The neutrinoless double beta (0u(3(3) decay, which can 
only occur by violating the conservation of the total lep- 
ton number, if observed, will unravel physics beyond 
Standard Model (SM), and will represent a major mile- 
stone in the study of the fundamental properties of neu- 
trinos Recent results from neutrino oscillation 
experiments have convincingly demonstrated that neu- 
trinos have mass and they can mix 0-0- The neutri- 
noless double beta (Ov/3/3) decay is the most sensitive 
process to determine the absolute scale of the neutrino 
masses, and the only one that can distinguish whether 
neutrino is a Dirac or a Majorana particle. A key ingre- 
dient for extracting the absolute neutrino masses from 
0v(3(3 decay experiments is a precise knowledge of the 
nuclear matrix elements (NME) for this process. Since 
most of the (3(3 decay emitters are open shell nuclei, many 
calculations of the NME have been performed within the 
pnQRPA approach and its extensions [To[-[2l|. However, 
the pnQRPA calculations are very sensitive to the vari- 
ation of the so called g pp parameter (the strength of the 
particle-particle interactions in the 1 + channel) [lQ|-[l2|, 
and this drawback still persists in spite of various im- 
provements brought by its extensions fTpl ]- p~8(] « including 
higher-order QRPA approaches [l^|-|2lj|. The outcome 
of these attempts was that the calculations became more 
stable against g pp variation, but at present there are still 
large differences between the values of the NME calcu- 
lated with different QRPA-based methods, which do not 
yet provide a reliable determination of the two-neutrino 
double beta (2u(3(3) decay half-life. Therefore, although 
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the QRPA methods do not seem to be suited to predict 
the 2v(3(3 decay half-lives, one can use the measured 2v(3(3 
decay half-lives to calibrate the g pp parameters that are 
further used to calculate the 0v(3(3 decay NME [22[. An- 
other method that recently provided NME for most 0v(3(3 
decay cases of interest is the Interactive Boson Model 
[23I ]. Given the novelty of these calculations, it remains 
to further validate their reliability by comparison with 
experimental data. 

On the other hand, the progress in computer power, 
numerical algorithms, and improved nucleon-nucleon ef- 
fective interactions, made possible large scale shell model 
calculations of the 2v(5(5 and 0v(3(3 decay NME 
The main advantage of the large scale shell model cal- 
culations is that they seem to be less dependent on the 
effective interaction used, as far as these interactions are 
consistent with the general spectroscopy of the nuclei in- 
volved in the decay. Their main drawback is the limi- 
tation imposed by the exploding shell model dimensions 
on the size of the valence space that can be used. The 
most important success of the large scale shell model cal- 
culations was the correct prediction of the 2v(3(3 decay 
half-life for 48 Ca 0, [27[. In addition, these calcula- 
tions did not have to adjust any additional parameters, 
i.e. given the effective interaction and the Gamow- Teller 
quenching factor extracted from the overall spectroscopy 
in the mass-region (including beta decay probabilities 
and charge-exchange form factors), one can reliably pre- 
dict the 2v(3(3 decay half-life of 48 Ca. 

Clearly, there is a need to further check and refine these 
calculations, and to provide more details on the analy- 
sis of the NME that can be validated by experiments. 
We have recently revisited [28[ the two neutrino double 
beta decay of 48 Ca using two recently proposed effective 
interactions for this mass region, GXPF1 and GXPF1A, 
and we explicitly analyzed the dependence of the double- 
Gamow- Teller sum entering the NME on the excitation 
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FIG. 1: (Color online) Comparison of the neutron and pro- 
ton occupation probabilities between GXPF1 interaction and 
FPD6 interaction. 



TABLE I: Neutron and proton occupation probabilities for 
the nuclei involved in the decay. 



Nucleus (N/P) 


Interaction 


O/7/2 


lP3/2 


O/5/2 


lPl/2 


48 CaN 


GXPF1 


7.883 


0.073 


0.033 


0.011 


48 Ti N 


GXPF1 


5.545 


0.237 


0.167 


0.051 


48 Ti p 


GXPF1 


1.846 


0.110 


0.033 


0.011 


4 *CaN 


GXPF1A 


7.892 


0.067 


0.032 


0.009 


48 Ti N 


GXPF1A 


5.535 


0.248 


0.168 


0.048 


48 Ti p 


GXPF1A 


1.839 


0.119 


0.032 


0.010 


48 CaN 


KB3 


7.800 


0.0706 


0.105 


0.024 


48 Ti N 


KB3 


5.422 


0.266 


0.248 


0.064 


48 Ti p 


KB3 


1.770 


0.120 


0.089 


0.022 


4 «CaN 


KB3G 


7.795 


0.070 


0.112 


0.024 


48 Ti N 


KB3G 


5.416 


0.260 


0.263 


0.061 


48 Ti p 


KB3G 


1.763 


0.120 


0.097 


0.021 


48 CaN 


FPD6 


7.693 


0.161 


0.117 


0.029 


48 Ti N 


FPD6 


5.253 


0.369 


0.310 


0.068 


48 Ti p 


FPD6 


1.673 


0.196 


0.101 


0.031 



the high-energy gamma and beta radiation emitted in 
this process could help eliminate most of the background 
noise. However, the small natural abundance of this iso- 
tope, 0.187%, increases the difficulty of an experimental 
investigation, although new, improved separation tech- 
niques were recently proposed [32]. In addition, previous 
calculations [26|, [33[ suggest that its NME is a factor 4-5 
smaller that those of other (3(3 emitters, such as 76 Ge and 
82 Se [46]. Since these calculations were reported, it was 
shown that the short range correlations might not have 
such a dramatic effect on the NME HHGH as previously 
thought, and it was also shown that higher order terms 
in the nucleon currents could be important [35[. In the 
present paper we take into account all these new devel- 
opments in the analysis of the NME for the Ov/3/3 decay 
of 48 Ca, and we study the NME dependence on the short 
range correlations parameters, on the finite size param- 
eter, on the average energy of the intermediate states, 
and on the effective interaction used for the many-body 
calculations. 



II. THE NEUTRINOLESS DOUBLE BETA 
DECAY MATRIX ELEMENT 

The 0v(3(3 decay (Z, A) (Z + 2, A) + 2e~ requires 
the neutrino to be a Majorana fermion, i.e. it is identical 
to the antineutrino. Considering only contributions from 
the light Majorana neutrinos @, the 0v(3(3 decay half-live 
is given by 
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Ov |2 
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(1) 



Here, G\ v is the phase space factor, which depends on 
the 0iy(3(3 decay energy, Qpp and the nuclear radius. The 
effective neutrino mass, (ra/3/3), is related to the neutrino 
mass eigenstates, ra&, via the lepton mixing matrix, U e k, 



energy of the 1 + states in the intermediate nucleus 48 Sc. 
This sum was recently investigated experimentally [291 ], 
and it was shown that indeed, the incoherent sum (us- 
ing only absolute values of the Gamow- Teller matrix ele- 
ments) would provide an incorrect NME, validating our 
prediction. We have also corrected by several orders of 
magnitude the probability of transition of the g.s. of 
48 Ca to the first excited 2 + state of 48 Ti. Future experi- 
ments on double beta decay of 48 Ca (CANDLES [U and 
CARVEL [31]) may reach the required sensitivity of mea- 
suring such transitions, and our results could be useful 
for planning these experiments. 

In the present paper we continue our investigations of 
the double beta decay of 48 Ca, analyzing the 0u(3(3 decay 
NME. 48 Ca has the largest Q^- value of 4.271 MeV (next 
largest is that of 150 Nd decay, 3.367 MeV), which could 
contribute to an increased decay probability. In addition, 



M=iE m ^i • ( 2 ) 

k 

The NME, M 0zy , is given by 

M 0u = M% V T - Mf - M^ v , (3) 

where M% V T , and M% v are the Gamow-Teller (GT), 
Fermi (F) and tensor (T) matrix elements, respectively. 
These matrix elements are defined as follows: 

Mr = ^(0+|r_ ro r_„O«j0+) , (4) 

m,n 

where < ^ nn are Ov/3/3 transition operators, 
a = (GT, F, T), | Of > is the g.s. of the parent 
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TABLE II: Parameters for the SRC parametrization of Eq. 
(CD- 



TABLE III: Different contributions to the NME for the 
GXPF1A interaction with (E) = 7.72 MeV 



SRC 


a 


b 


c 


SRC 






Mt u 




Miller-Spencer 


1.10 


0.68 


1.00 


None 


0.556 


-0.219 


-0.015 


0.711 


CD-Bonn 


1.52 


1.88 


0.46 


Miller-Spencer 


0.465 


-0.141 


-0.014 


0.570 


AV18 


1.59 


1.45 


0.92 


CD-Bonn 


0.688 


-0.222 


-0.014 


0.845 










AV18 


0.634 


-0.204 


-0.014 


0.779 




No SRC 
■ M-S SRC 

CD-Bonn SRC 
AV18 SRC 



GXPF1 GXPF1A KB3 KB3G FPD6 

FIG. 2: (Color online) The dependence of the NME on the ef- 
fective interaction used and the short range correlation (SRC) 
model. M-S stands for Miller-Spencer. 

nucleus (in our case 48 Ca), and | 0^ > is the g.s. of the 
grand daughter nucleus (in our case 48 Ti). 

Due to the two-body nature of the transition oper- 
ator, the matrix element can be reduced to a sum of 
products of two-body transition densities TBTD and 
antisymmetrized two-body matrix elements, 

AC = J2 TBTD(j p j p ,,j n j n ,;J n ) 

jpjp' JnJn' Jir 

(j p j p ,;J"T | r_ 1 r_ 2 1 « 2 | j n jw;J 7r T) a , (5) 
where Of 2 are given by 

Cf 2 T = <7i • a 2 H GT {r) , 
0[ 2 = H F {r) , 
0\2 = I 3 (<?i ' ?) (°i -r) -ai- a 2 ] H T (r) . 



(6) 



The matrix elements of 0" 2 f° r ^ ne jj-coupling scheme 
consistent with the conventions used by modern shell 
model effective interactions are described in the Ap- 
pendix. 

To calculate the two-body matrix elements in Eq. (5) 
one needs the neutrino potentials entering into the radial 
matrix element (nl | H a | nl') in Eq. ([T3]) below. Follow- 
ing Ref. [35] and using closure approximation one gets, 



H a (r) 



2R r 

K Jo 



1 fa(Qr)^QrG a (q 2 )qdq , (7) 



where f F ,GT(qr) = jo(qr) and f T (qr) = 32{qr) are spher- 
ical Bessel functions, (E) is the average energy of the 



virtual intermediate states used in the closure approx- 
imation, and the form factors h a (q 2 ) that include the 
higher order terms in the nucleon currents are 



h F (q 2 ) 

h.GT(q 2 ) 



gW) 
9W) 



g\ 
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3 q 2 + m 2 3 \q 



h T (q 2 



2 g 2 M (q 2 ) g 2 
3 g\ 4m 2 p ' 



9 2 a(Q 2 ) 



9a 



1 



3 q 2 + m 2 3 \ q 2 

1 9 2 m(Q 2 ) I 2 
3 9a 4m2 ' 



(8) 



The gv,A,M form factors in Eq.® can include nucleon 
finite size effects, which in the dipole approximation are 
given by 



9v(q 2 ) = 

M 2 ) = (M 
9A(q 2 ) = 



gv 



2 ' 



(l + 4 2 /A 2 v) 
9M(q 2 ) - Hn)gv(q 2 ), 

^\ - 9A 

(l + ^ 2 /Ai) 2 ' 



(9) 



Here g v = 1, g A = 1.25, (/i p -/i n ) 3.7, Ay 850 MeV, 
and A A = 1086 MeV. 

The short range correlations are included via the corre- 
lation function f(r) that modifies the relative wavefunc- 
tions at short distances, 



i>nl(r)^[l + f(r)]i>nl(r) , 
where /(r) can be parametrized as |35l |. 

f(r) = -ce- ar2 (1 - br 2 ) . 



(10) 



(11) 



The radial matrix elements of H a between relative 
harmonic oscillator wavefunctions ijj n i(r) and ^ n ///(r), 
(nl | H a (r) | n'V), become 



/ r 2 dr^ nl {r)H a {r)^ v {r) [1 + /(r)] 2 (12) 
J o 



4 




<E> (MeV) 



FIG. 3: (Color online) The dependence of the NME on the av- 
erage energy of the intermediate states, (E), for the GXPF1A 
interaction. 



FIG. 4: (Color online) The dependence of the NME on the 
axial cutoff parameter A a for the GXPF1A interaction. 



Although the neutrino potentials are quite close to a 
Coulomb potential, the integrands in Eq. (7) are strongly 
oscillating, and the integrals require special numerical 
treatment. 

Having calculated the two body matrix elements, we 
developed a new shell model approach for computing the 
many-body matrix elements for Oupp transition, Eq. (5). 
This approach is briefly described in the Appendix. 



III. RESULTS 

In this study of the decay NME we used five dif- 

ferent effective interactions available for the shell model 
description of the p/-shell nuclei: GXPFjJl, GXPF1A 
[H, KB3 [H, KB3G 0, and FPD6 [U These ef- 
fective interaction were constructed starting from a G- 
matrix [42[ in the pf shell, which was further adjusted 
to describe some specific (but different) sets of experi- 
mental energy levels of some p/-shell nuclei. Although 
their elements are quite different, their predictions of the 
spectroscopic observables around A=48 are not very far 
apart. Recent experimental investigations [ijjl 53] of the 
nucleon occupation probabilities in 76 Ge and 76 Se and 
the subsequent theoretical analysis [35|, [37| highlighted 
the relevance of these observables to obtain an accurate 
description of the nuclear structure of the nuclei involved 
in the double beta decay. Fig. [1] compares the neutron 
and proton occupation probabilities in 48 Ca and 48 Ti for 
two different effective interactions, GXPF1 and FPD6. 
One can see very small differences between the results of 
the two interactions. One can come to the same conclu- 
sion when comparing the similar occupations probabili- 
ties for all five interactions given in Table H 

In the present calculations we considered both short 
range correlation (SRC) effects and finite size (FS) ef- 
fects. Although the radial dependence of the neutrino po- 



tential is very close to that of a Coulomb potential, many 
previous calculations [25|, [26|, [33|, [46[ took into account 
the SRC missing in the two-body product wave func- 
tions, via the Jastrow-like parametrization described in 
Eqs. (10)-(12). Until recently, the parameters a, 6, c used 
were those proposed by Miller and Spencer [45[, which 
have the effect of decreasing the NME by about 30%. 
Recently, [35] the SRC effects were revisited, using mod- 
ern nucleon-nucleon interactions, such as CD-Bonn and 
AVI 8, and it was found that the decrease of the relative 
wave functions at short distances is compensated by a rel- 
ative increase at larger distances, and the overall NME 
do not change very much when compared with the NME 
without SRC effects. Ref. [35j proposed a parametriza- 
tion of these results in terms of similar Jastrow-like cor- 
relation functions as in Eq. (10)- (11), the corresponding 
parameters being listed in Table HI1 In addition, Ref. [34[ 
introduced an effective Ov/3/3 operator that takes into ac- 
count the SRC effects and the contribution of the missing 
shells from the valence space using the general theory of 
effective interactions [421 ]. and found that the NME for 
the Ov/3/3 decay of 82 Se did not change significantly when 
compared with the result of the "bare" operator. 

Fig. [2] shows our NME for all five effective interac- 
tions, for all three SRC set of parameters given in Table 
HI! as well as for no SRC. One can see that the semi- 
quantitative analysis above is reflected in the dependence 
of NME on the choice of SRC. The results do not show 
significant dependence on the effective interaction used, 
although one can see a 20% spread of NME for the same 
choice of SRC. All NME reported here contain the higher 
order terms described in Eqs. (7)- (9). A comparison 
with the NME calculated without the higher order terms 
in the potential will be reported elsewhere. To be consis- 
tent [43 with the calculation of the phase factor, G^, we 
used R = 1.2 A 1 / 3 fm in Eq. (7). Our choice for hco pa- 
rameter entering the harmonic oscillator wave functions 
was 45A -1 / 3 — 25A -2//3 , which was shown to provide bet- 
ter shell model description of observables than the simple 
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41A _1/3 anstaz. Table Hill shows the GT, F, and T con- 
tributions to the overall NME for all SRC choices, when 
the GXPF1A interaction was used. One can see that all 
theses contributions add coherently in Eq. (3), and that 
the tensor contribution is negligible in all cases. 

Fig. [3] shows the dependence of the NME of the average 
energy of the intermediate states. Varying (E) from 2.5 
MeV to 12.5 MeV one gets less than 5% variation in the 
NME. This result suggests that closure approximation is 
quite good, although a direct study might be necessary to 
find out the exact magnitude of the error. All other NME 
results reported here used (E) = 7.72 MeV. Fig. 0] shows 
the dependence of NME of the axial cutoff parameter 
A a, while Ay is kept fixed at 850 MeV. The variation 
is within 5%, indicating a weak dependence of the FS 
cutoff parameters. The FS effects were implemented via 
the cutoff parameters Ay and A^ in the form factors 
given in Eq. (9). In most of our results, except those 
in Fig. [4]we use the same FS cutoff parameters as in 
Refs. IfllHl, i.e. Ay = 850 MeV and A A = 1086 MeV. 
In Fig. 2]we present the dependence of the NME on A^, 
while Ay is kept at its nominal value. As with the (E) 
dependence, the results varies by less than 5%. 

IV. CONCLUSIONS AND OUTLOOK 

In conclusion, we presented a new shell model anal- 
ysis of the nuclear matrix elements for the neutrinoless 
double beta decay of 48 Ca. We included in the calcu- 
lations the recently proposed higher order terms of the 
nucleon currents, three old and recent parametrization 
of the SRC effect, FS effects, intermediate states en- 
ergy effects, and we treated careful few other parame- 
ters entering the into the calculations. We found very 
small variation of the NME with the average energy of 
the intermediate states, and FS cutoff parameters, and 
moderate variation vs the effective interaction and SRC 
parametrization. Our overall average NME using all val- 
ues presented in Fig. [2] is 0.86, although the elimina- 



tion of the Miller-Spencer SRC parametrization would 
significantly increase this value. We estimate the error 
due to the effects studied here to be about 18%. Using 
the present value of the NME and the recommended J47[ 
phase-space factor G\ v — 6.5 x 10 -14 yr~ x one can con- 
clude that a future measurement of the Ov/3/3 decay half- 
life of 10 26 yr, which seems to be the limit imposed by 
the present energy resolution of the CANDLES detector 
[32( | (see also Fig. 21 of Ref. could detect a neutrino 

mass (mpp) of about 230 ± 45 meV. New improvements 
in the detector technology could further reduce this limit. 

We believe that our analysis covered the most impor- 
tant effects relevant to the accuracy of the NME for the 
double beta decay of 48 Ca. The successful prediction of 
the 2vj3j3 decay half-live of 48 Ca using the shell model ap- 
proach in the same model space suggests that the Oz^/3/3 
decay half-life could be reliably predicted. Our analysis 
indicates that closure approximation is accurate at the 
level of 5% error. However, a direct comparison with the 
sum on intermediate states, calculated within the shell 
model, would be more reassuring. 

As in all other Ov/3/3 decay calculations reported so 
far, no quenching factors have been used. This is prob- 
ably the least investigated potential effect on the OvP/3 
decay NME. It is possible that the GT-like operators be 
quenched even if they have finite range. No definitive 
conclusion exists about these effects, although Ref. [34| 
seems to indicate that they might not play a major role. 
One could try getting some insight into this problem by 
investigating the fictitious 0^/3(3 decay of a light nucleus, 
such as 16 Be, and increasing the valence space one can 
study the changes of different contributions to the NME. 

V. APPENDIX 

The matrix element of Of 2 f° r non-antisymmetrized jj- 
coupling two-particle states can be decomposed into 
products of reduced matrix elements of the spin, relative 
and center of mass operators, 



(ni/iji,n 2 Z 2 j 2 ; J*T = 1 | 0? 2 | n[l[j[,r^l' 2 j' 2 ;J"T = l) = 

x {nl,NL\n 1 l 1 ,n 2 l2) x (n'l,NL\n[l[,n 2 r 2 ) x {nl\ H a (r)\n'l). (13) 

nn'INL 

where 



h\{h),h\{h)\\\{S)Mh{X) 



V(2ji + l)(2j 2 + l)(2S+l)(2A 




(14) 



and the last factor is a 9j symbol. As mentioned in the tent with the convention considered into the derivation 
text, the order of the spin-orbit coupling must be consis- 
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of the effective interaction used for the many body cal- 
culations (I + s in most cases). The reduced spin matrix 
elements are 

(^ll^i-^H^s) = \2S{S + l) -3] , 

(^IWI^) = VWTT. (15) 

The expressions for the matrix elements of the ten- 
sor operator will be given elsewhere. Our calculations 
show that the tensor term gives a small contribution to 
the NME in most cases. The higher order terms in the 
nucleon currents, however, decrease the overall NME by 
about 20-25%. The antisymmetrized form of the two- 
body matrix elements can be obtained using 

{3p3p>;J*T | *_if_ 2 Of 2 | U3n>\ J^T) a = 
1 

l(j P jp>;J n T | *_it- 2 0? 2 | j n j n >]J*T) - 
(_l)in+i n ,+J ( jpjp ,.j* T | *_i*_ 2 0? 2 | j n ,j n] J*T>] (16) 

Having the two-body matrix elements ready, one can 
calculate the NME in Eq. ([5]) if two body transition den- 
sities TBTD (jpjp> , jnjri', Jtt) are known. Most of the 
shell model codes do not provide two body transition 
densities. One alternative approach is to take advantage 
of the isospin symmetry that most of the effective in- 
teractions have, which creates wave functions with good 
isospin. The approach described below works also when 
the proton and neutron are in different shells. If the 
above conditions are satisfied, one can transform the two 
body matrix elements of a AT = 2 operator using the 
Wigner-Eckart theorem, from AT Z = -2 to AT Z = 0, 
which can be further used (see below) to describe transi- 
tions between states of the same nucleus. Denoting 

<o£?r= 2 - 2 > = (it) 

(T = 1 T z = -1 | 0££= = 2 _ 2 | T = 1 T z = 1) , (18) 
one gets for the AT Z — the two body matrix elements 

< PP I O^7 = 2 I PP >= 
< nn | Oatw) I nn >= 
<0&=?_ 2 >xCiil/ClJ_\, (19) 

and 



< pn T = 1 | O££L 2 \pnT=l >= 

< OaT z = -2 > X Q) o/^l -2 -1j (20) 

where C\ z . \^ T are isospin Clebsch-Gordan coeffi- 
cients. The transformed matrix elements in Eqs. ([T9|) 
and (|20|) preserve spherical symmetry and they can be 
used as a piece of a Hamiltonian, H** , that violates 
isospin symmetry. 

One can then lower by 2 units the isospin projection of 
the g.s. of the parent nucleus that has the higher isospin 
T>, 48 Ca in our case, thus becoming an isobar analog 
state of the grand daughter nucleus that has isospin T< = 
T> — 2, 48 Ti in our case. Denoting by | 0^ T> > the 
transformed state, one can now calculate the many body 
matrix elements of the transformed 0v[3[3 operator, 

Ml%T z = r<) = (o+ T< I Hpp | 0+ T>) . (21) 

Choosing | 0^ T> > as a starting Lanczos vector, and 
performing one Lanczos iteration with H%p one gets 

Hfa | 0+ T> >= ax | 0+ T> > +b x \ L t >, (22) 

where | L\ > is the new Lanczos vector. Then, one can 
calculate the matrix elements in Eq. (f2Tj) 

M°J{T Z = T<) = 6i (0+ T< | Li) . (23) 

The transition matrix elements in Eq. ([5]) can then be 
recovered using again the Wigner-Eckart theorem, 

AC = M° a »(T z =T < )x C T T l _ 2 2 t < k IC T T > I T T < . (24) 

Although it looks complicated, this procedure is rather 
easy to implement. The transformation of the g.s. of the 
parent to an analog state of the grand daughter can be 
performed very quickly, and one Lanczos iteration repre- 
sents a small load as compared with the calculation of the 
g.s. of the grand daughter. The additional calculations 
described in Eqs. (22)- (24) require smaller resources than 
those necessary to calculate the TBTDs. 
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